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Radio pulsars with millisecond spin periods are thought to have been spun up 
by transfer of matter and angular momentum from a low-mass companion 
star during an X-ray-emitting phase. The spin periods of the neutron stars in 
several such low-mass X-ray binary (LMXB) systems have been shown to be 
in the millisecond regime, but no radio pulsations have been detected. Here we 
report on detection and follow-up observations of a nearby radio millisecond 
pulsar (MSP) in a circular binary orbit with an optically identified companion 
star. Optical observations indicate that an accretion disk was present in this 
system within the last decade. Our optical data show no evidence that one 
exists today, suggesting that the radio MSP has turned on after a recent LMXB 
phase. 

The fastest-spinning radio MSPs are thought to be formed in systems containing a neutron 
star (NS) and a low-mass (< 1M ) companion star fj). Mass transfer occurs when matter over- 
flows the companion's Roche lobe (T21), forms an accretion disk around the NS, and eventually 
falls onto its surface, producing bright X-ray emission ([/][3]). 

Radio emission that would otherwise be produced by a rapidly rotating magnetic NS is 
thought to be quenched during active accretion by the presence of ionized material within the 
pulsar's light cylinder. However, for a sufficiently low accretion rate, infalling material may 
be halted outside the light cylinder by magnetic pressure, presumably allowing the MSP's radio 
emission to turn on. At this point, the MSP's electromagnetic emission and particle wind should 
irradiate the disk and companion, driving mass out of the system. This activation as a radio 
pulsar, a possible end to the accretion phase, is poorly understood, in part because no radio 
MSP has yet been observed to turn on in such a system, although there is indirect evidence that 
this may occur (|?][5][6]). 

The radio source FIRST J102347.67+003841.2 (hereafter J1023) appears optically as a V ~ 
17.5 magnitude star with a mid-G (solar-type) spectrum and mild 0.198-day orbital variability 
0. However, in observations from May 2000 to Dec. 2001, J1023 had a blue spectrum with 
prominent emission lines (absent in quiescence) and exhibited rapid flickering by ~ 1 
magnitude (ph. This optical behavior is typical of an accretion flow, and the double-peaked 
nature of the lines (|9|) suggests an accretion disk specifically. This history led to the system's 
classification as an accreting white dwarf §8§ or NS Q binary. For a more detailed history 
see the Supporting Online Material (SOM). Since 2002 the system has shown no further sign 
of accretion, and our optical spectroscopic observations of Dec. 23 and 25, 2008 (see SOM) 
confirm the continued absence of emission lines. Thus J1023 is currently in a quiescent state. 

We found PSR J1023+0038, a bright MSP with a spin period of 1.69 ms at a dispersion 
measure (DM) of 14.33 pc cm -3 , as part of a 350-MHz pulsar survey carried out with the 
Robert C. Byrd Green Bank Telescope (GBT) in 2007 (see SOM for details). 

Phase-coherent timing of the radio pulsations at the GBT, the Arecibo Observatory, and 
Parkes Observatory has allowed high-precision measurements of the binary parameters (see 
Tableland SOM for details). Together with the much smaller position uncertainties implied 
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by the telescope beam sizes of the timing observations, these binary parameters unambiguously 
place this pulsar in the J1023 system, confirming the prediction of (7) that J1023 has a NS 
primary. 

Combining these parameters with optical radial velocity measurements §7\, we find that the 
pulsar is a factor of 7.1 ± 0.1 more massive than its companion. The determination of the 
individual masses depends only on the inclination angle i of the orbital plane from the plane of 
the sky. The pulsar's Keplerian mass function (m c sin i) 3 / [m p + m c ) 2 is 1.1 x 10~ 3 M , where 
m c and m p are the masses of the companion and the pulsar, respectively. For a NS mass in the 
range ~ 1.0 to ~ 3.0 M ( [TO] ), i is restricted to the range ~ 53° to ~ 34°, and the companion 
mass is restricted to ~ 0.14 to ~ 0.42M Q . 

Stellar models, based on the assumption (reasonable in light of the evidence for a recent 
accretion flow) that the companion is filling its Roche lobe, show that it is possible for the com- 
panion star to exhibit the observed orbital modulation of color, magnitude, and radial velocity, 
if the primary has an isotropic luminosity of ~ 2 L (7). These models also imply i < 55°, 
which is compatible with the range obtained from theoretical NS mass limits. Similar models 
account for the light curve modulation during quiescence of the low-mass companion to SAX 
Jl 808.4—3658 §6§, an accreting X-ray pulsar in a 2 hr orbit (e.g. [TT][T2]), and which has been 
suspected of harboring an MSP when in quiescence ([? 13 5]). 

The brightness, radial velocities, and presumed Roche lobe filling companion of J1023 im- 
ply that the distance to the system is 0.9 kpc/ sini; for the allowed range of i, the distance is 
between ~ 1.1 and ~ 1.6 kpc. A pulsar mass of 1.4M implies a distance of 1.3 kpc. For 
comparison, given the observed DM of 14.325 pc cm -3 the standard Galactic free-electron 
density model ( |7?| ) predicts a distance of 0.6 kpc, consistent with our range given its known 
uncertainties. 

Ordinarily one can estimate a pulsar's spin-down luminosity from its observed spin-down 
P, but for MSPs this is subject to contamination by several factors, including accelerations of 
J 1023 and the solar system barycenter in the Galactic gravitational potential. From the measured 
optical proper motion /i of J 1023 (Table [TJ we expect a significant positive contribution to 
P = Pdfi 2 /c ~ 1.8 x 10 -21 (<i/1.3 kpc) due to the Shklovskii effect §L5§. While in some binary 
systems the accelerations can be constrained independently using orbital period variations, in 
J 1023 we observe very large orbital period variability (binary period derivative ~ 3 x 10~ 10 ), 
almost certainly due to classical tidal torquing from gas motion in the extended envelope of the 
companion star. The near- zero eccentricity of the pulsar's orbit (Table [T]) also suggests that the 
system has undergone tidal circularization. 

This picture of a Roche-lobe- filling companion is supported by evidence for the presence 
of ionized material in the system. Our broad-band radio observations reveal substantial DM 
variations at certain orbital phases (Fig. [T] G) as well as smaller DM variations on time scales 
of minutes throughout the orbit. Moreover, we observe regular eclipses ranging from very brief 
at 3000 MHz to most of the orbit in our 150 MHz observations with the Westerbork Synthesis 
Radio Telescope (Fig. [T] A-F). We also see very brief eclipses at all orbital phases (Fig. [T]E 
in particular). Given the established range of i, the line-of-sight between the pulsar and the 
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Earth will not intersect the Roche lobe of the companion at any point in the orbit. If we assume 
that during the DM variations near orbital phase 0.4 (Fig. [T]ii and iii) the obstructing plasma is 
uniform, has a plasma frequency of 700 MHz (based on the loss of the 700 MHz signal at this 
phase), and contributes the 0.15 pc cm~ 3 excess DM, we predict a thickness of 3 x 10 4 km (thin 
compared to the companion size of ~ 3 x 10 5 km). This suggests that a thin but dense layer of 
material, perhaps a shock due to the pulsar wind meeting either winds from the companion or 
material overflowing the companion's Roche lobe, crosses the line of sight at this orbital phase. 

Similar eclipses and DM variations have been seen in the Galactic "black widow" systems, 
tight NS binary systems in the Galactic field with companions an order of magnitude less mas- 
sive than that of J1023 (e.g. [16]). Such eclipses have also been seen in several radio MSPs in 
globular clusters (e.g. [T7][T8| ). In both cases they have been attributed to the presence of gas 
flowing out from the irradiated companion, as is probably the case for J 1023. However J 1023 
differs from all other known radio MSPs in that there is evidence for a recent accretion disk. 

For an accreting NS the accretion luminosity is proportional to the accretion rate, which is 
unknown for J 1023. But if the matter is to reach the NS surface, it must overcome magnetic 
pressure at the Keplerian corotation radius, so there is a minimum accretion rate and hence a 
minimum luminosity for a given NS magnetic field (19). For a magnetic field equal to our 
timing-derived upper limit (Table [I]) this minimum luminosity for J1023 is 1.1 x 10 37 erg s _1 
(using formulae from [19]). Such a high luminosity can be clearly ruled out by our analysis of 
archival Rossi X-ray Timing Explorer all-sky monitor data (see SOM for details): the average 
annual X-ray luminosity for each year 1996-2008 was less than 4.8 x 10 33 erg s -1 (<i/1.3 kpc) 2 . 
The variable nature of LMXB X-ray emission cannot explain this low luminosity: on Feb 1, 
2001, when an optical emission line spectrum was observed (9), the average flux was less than 
2.4 x 10 34 erg s _1 (<i/1.3 kpc) 2 . If accretion occurred during the active phase, the annual upper 
limit would imply that the magnetic field of the NS in J1023 was less than 6 x 10 6 G, smaller 
than that of any known MSR 

It therefore seems more likely that infalling matter did not reach the NS surface, but instead 
underwent what is known as "propeller- mode accretion" (20): infalling matter entered the light 
cylinder but was stopped by magnetic pressure outside the corotation radius, which prevented 
it from falling further inward. This process also has a minimum accretion rate and minimum 
luminosity: if the infalling matter does not reach the light cylinder, the radio pulsar mecha- 
nism will presumably become active. If this occurs no stable balance exists between outward 
radiation and wind pressure and ram pressure of infalling material, and the disk will be cleared 
from the system (]79j). The minimum luminosity for propeller-mode accretion is substantially 
lower than that for standard accretion: our timing-derived upper limit on the magnetic field 
(Table [I]) implies a minimum luminosity of just 2.7 x 10 32 erg s _1 (using formulae from 19). 
Thus it seems likely that during J1023's active phase, the mass transfer rate was high enough 
for propeller-mode accretion, but not high enough for material to reach the NS surface. 

In this scenario, a small drop in the mass transfer rate would clear the disk from the system 
and return J 1023 to its current quiescent state. Overflowing matter from the companion would 
then encounter a shock near the inner Lagrange point and be carried out of the system ([79|), 
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possibly explaining the observed variable hard power-law X-ray spectrum of J1023 (21) and 
the presence of gas in the system. Such shocks have been suggested to explain the variable 
power-law X-ray emission in 47 Tuc W (|22) and PSR J 1740-5340 in NGC 6397 ((25]), which 
also have low-mass companions, and have previously been compared with SAX Jl 808.4— 3658, 
though with no evidence for the recent presence of an accretion disk. 

Despite the apparent resemblance of J1023 to 47 Tuc W and PSR J1740-5340, the latter 
two systems (and indeed all other currently known eclipsing pulsars with several-tenths M 
companions) reside in globular clusters, and are likely to have acquired their current companions 
in exchange interactions after the pulsars were "recycled" d2?l). J1023 is the only known highly 
recycled (the fifth fastest known) MSP in the field of the Galaxy with both a non-degenerate 
companion and an orbit that has been circularized through tidal interactions. A globular cluster 
origin for J 1023 is extremely unlikely due to its large distance from the nearest globular cluster 
as well as from the Galactic bulge. The evidence points to J 1023 's having been recycled by its 
current companion, which has not yet completed the transformation to a white dwarf. 

The observed transition of J1023 suggests that it is in a bistable state: for certain rates 
of Roche lobe overflow, if the radio pulsar mechanism is quenched, propeller-mode accretion 
can occur, but if the radio pulsar mechanism is active, that same mass accretion rate cannot 
overcome the radiation pressure and no accretion occurs. Should the mass transfer rate of 
J 1023 rise sufficiently, then, it may enter another LMXB phase: a disk will form, the radio 
emission may be quenched, and the X-ray luminosity may increase dramatically, due to either 
propeller-mode accretion, with a net spin-down of the pulsar, or even brighter accretion onto 
the surface, with a net spin-up. 
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Parameter 


Value 


Right ascension (a; J2000) 


10 h 23 m 47.687(3) s 


Declination (5; J2000) 


00°38'41.15(7)" 


Proper motion in a (fx a ) 


10(1) mas yr _1 


Proper motion in 5 


— 16(2) mas yr _1 


Epoch (MJD) 


54802 



Dispersion Measure 

Pulsar period (P) 

Pulsar period derivative (P) 

Orbital period (P b ) 

Orbital period derivative (P&) 

Orbital period second derivative (P&) 

Time of ascending node (MJD) 

asini 

Eccentricity 



14.325(10) pc cm" 3 
1.6879874440059(4) ms 
1.2(8) x 10~ 20 
0.1980962019(6) d 
2.5(4) x 10" 10 
-5.21(14) x 10~ n s- 1 
54801.97065348(9) 
0.3433494(3) lt-s 
< 2 x 10~ 5 



1600 MHz flux density 
Spectral index (a; S oc v~ a ) 
Surface magnetic field 
Spin-down luminosity 



~ 14 mJy 

~ -2.8 

< 3 x 10 8 G 

< 3 x 10 35 



erg s 



Table 1 Pulsar parameters. Parameters listed in the top section are held fixed in the timing fit 
that produces the parameters listed in the middle section (SOM). The position and proper mo- 
tion values are from the USNO NOMAD optical catalog ([25]). With the exceptions of DM, P, 
and eccentricity, uncertainties on timing quantities are twice the formal la errors returned by 
the pulsar timing analysis program. The DM of the source varies substantially: around orbital 
phase 0.3 it occasionally increases by ~ 0.15 pc cm -3 , and we observe apparent orbit-to-orbit 
variations as large as ~ 0.01 pc cm~ 3 . Our estimate of the eccentricity is approximate because 
DM variations limit the orbital coverage of usable timing data. Our estimate for P is approxi- 
mate because it is highly covariant with pulsar position; it was estimated by fitting for P while 
varying the position according to the optical uncertainties listed above. We have not subtracted 
Shklovskii or galactic accelerations. The flux density and spectral index are estimated using 
standard system temperature and gain values for the GBT and the Parkes telescope. The flux 
density is subject to substantial variability intrinsic to the pulsar (we have observed brighten- 
ing by a factor ~ 4 within minutes), as well as interstellar scintillation; for this reason we 
estimated the spectral index using the simultaneous 600 MHz and 3000 MHz Parkes observa- 
tions. Spin-down luminosity and magnetic field upper limits are calculated from the 3cr upper 
limit on P using the standard formulae E = 3.95 x 10 31 erg s _1 (P/10~ 15 )(P/l s)~ 3 and 

B = 3.2 x 10 19 G \fp~P. 
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Figure 1 Frequency dependence of eclipses, DM variations, and pulse profiles. (A-F) show 
flux density as a function of pulse phase and orbital phase at 3000, 2000, 1600, 700, 350, and 
156 MHz, respectively. (G) shows DM, as estimated from the 1600 MHz observations (see 
SOM for details). Orbital phase is defined to be zero at the pulsar's ascending node, so that the 
companion passes closest to our line of sight at orbital phase 0.25. To the right of (A) through 
(F) are pulse profiles at the respective bands; the instrumental smearing time is indicated with a 
horizontal bar in each panel. The exception is that the profile (C) is instead based on a 1410 MHz 
observation with higher time resolution. Note that the eclipse is nearly absent at 3000 MHz, but 
is longer at 700 MHz than at 1600 MHz. At low frequencies, random short eclipses, indicated 
by (iv) and (v), are also visible. Pulse phase is as predicted by a phase-coherent timing solution, 
so that large short-term timing variations are visible as vertical motions of the pulse peak. Note 
in particular the pulse arrival time variations at eclipse ingress (i) and egress (ii); the variation 
at egress appears to be due to the substantial DM variation (iii) visible in (G). 
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Supporting Online Material 



History of observations of J1023: The source FIRST J102347.67+003841 was detected as 



a point source in 1998 August in the FIRST VLA 1.4 GHz continuum radio sky survey (26). 
Three observations over the course of a week showed that the flux density of J1023 varied 
from less than a few mJy to roughly 7 mJy §8§. The source was also identified with a V ~ 17.5 
magnitude star having a history of observations back to a 1952 Palomar Observatory Sky Survey 
observation ( |25] ). In 1999 March, a Sloan Digital Sky Survey filter-photometry observation 
showed that this star had colors consistent with those of a G star (|9]). 

However, observations at optical wavelengths showed a dramatic change: spectra taken in 
May 2000 @ were substantially bluer, and had prominent emission lines, leading to classifica- 
tion of J 1023 as a cataclysmic variable. Further spectroscopic observations (|9][#|) confirmed this 
spectrum, and revealed emission lines with double peaks and orbital radial velocity modulation 
as late as December 2001. Optical photometry ft8h showed flickering, varying by ~ 1 mag on 
timescales from ten seconds to one day. 

Optical observations after December 2001 showed a return to the previous state: a May 2002 
observation ( [27] ) showed a G-like spectrum without emission lines, and photometry taken in late 
2002 and early 2003 d27|) shows a smooth light curve with a modulation period of 0.198 days, 
having amplitude ~0.3 mag in V, though with some low-level flickering. Observations in 2004 
found a smooth light curve consistent with that of 2003, though with no flickering, and 
a spectrum matching that of a mid-G star. Radial velocity measurements showed a regular 
Doppler shift with 0.198 day period and peak velocity 268 ± 4 km s _1 §7§. Light curve and 
color-variation modelling indicated that the primary was too massive and too luminous to be a 
white dwarf (given the observed red optical spectrum), so it was suggested that the source was 
probably an LMXB with a NS primary, now in quiescence (TTj). This idea was supported by 
XMM-Newton observations in 2004, which found a hard power-law X-ray spectrum ([27]). 

Further details of all these observations are compiled in Table [2] 

Our observations of J1023: In mid-2007, we carried out a 350-MHz pulsar survey with 
the Robert C. Byrd Green Bank Telescope (GBT) in West Virginia. This survey covered ~ 
12, 000 square degrees of sky in the declination range —21° to +26°, and used the Spigot pulsar 



autocorrelation spectrometer (28) to synthesize a filterbank sampling a 50 MHz band with 2048 
frequency channels every 81.92 fxs. Details of the survey and ongoing data analysis procedures 
will be presented elsewhere. In the data from June 28, 2007, we found PSR J 1023+0038, a 
bright MSP (mean flux density ~ 75 mJy in the 140 s discovery observation) with a spin period 
of 1.69 ms at a dispersion measure (DM) of 14.33 pc cm -3 . 

We obtained follow-up radio observations; all these observations are summarized in Ta- 
ble^ These observations used several different instruments and backends, in several different 
radio-frequency bands. Where possible, we observed simultaneously with both coherent dedis- 
persion backends with online folding (for precise timing) and incoherent dedispersion backends 
with offline folding (for larger bandwidths and to retain single-pulse and transient information). 
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For these observations we used the GBT, the Arecibo Observatory (AO) in Puerto Rico, the 
Westerbork Synthesis Radio Telescope (WSRT) in the Netherlands, and the Parkes telescope 
in Australia. The GBT Pulsar Spigot is an autocorrelation spectrometer (pSb. The Green Bank 
Ultimate Pulsar Processing Instrument (GUPPI) is a new digital filterbank for the GBT The 
Wideband Arecibo Pulsar Processors (WAPPs) are a set of autocorrelation spectrometers usable 
on independent frequency bands (29). The Parkes 10/50 cm dual-band receiver feeds a set of 1- 
bit analog filterbanks. The Arecibo Signal Processor (ASP) and the Green Bank Astronomical 
Signal Processor (GASP) are systems that use a computer cluster to coherently dedisperse ( |J0| ) 
and fold incoming data using the best available ephemeris. Further details of the operation of 
ASP and GASP and the reduction of their output data can be found in (31 ) and (32). We used 
the Pulsar Machine II (PuMa-II) coherent dedispersion backend at the WSRT (an interferom- 
eter we operated in phased-array mode). The precise modes we used for each observation are 
described in the footnotes to Table [3j 

In addition to the radio observations described above, we took an optical spectroscopic 
observation Dec. 23 and 25, 2008 with the Goodman high-throughput spectrograph on the 
Southern Astrophysical Research telescope. Since no flux calibration observation was taken, 
we produced a spectrum by averaging our seven 900-s exposures and dividing by the best-fit 
continuum. The resulting spectrum (Figure|2]) confirms the continued absence of emission lines, 
and hence quiescent state of the system. 



X-ray upper limit: In order to improve on the upper bound of (21 ) on J1023's X-ray bright- 
ness during its active phase, we searched the archival Rossi X-ray Timing Explorer All-Sky 
Monitor data for the known position of J1023. This data covers the energy range 2-10 keV for 
all of the span 1996-2006 with the exception of from brief periods when the source passed close 
to the Sun. The average annual flux for each year 1996-2008 must be beneath the systematic 
upper limit of 2.4 x 1CT 11 erg s _1 cm~ 2 (4.8 x 10 33 erg s -1 (<i/1.3 kpc) 2 ). Moreover, since 
LMXBs are variable sources, we searched individual days as well, obtaining varying upper 
limits; in particular, the average X-ray flux during Feb 1, 2001, when an optical emission-line 
spectrum was observed (9), did not exceed 1.2 x 1CT 10 erg s _1 cm -2 (corresponding to a 
luminosity of 2.4 x 10 34 erg s -1 (d/1.3 kpc) 2 ). 



Details of Main Text Figure [T| The data shown in (A) and (D) were taken simultaneously 
with the Parkes 10/50 cm dual-band receiver and analog filterbank on Nov. 24, 2008; the data 
shown in (B) were taken with the GBT Pulsar Spigot on Oct. 18, 2008; the data shown in (C) 
and (G) are the same 1600 MHz observation taken with the GUPPI and the GBT Pulsar Spigot 
respectively on Nov. 1, 2008; and the data in (E) and (F) were taken with the Westerbork 
Synthesis Radio Telescope using the PuMa-II coherent dedispersion backend on Nov. 11, 2008 
and Dec. 31, 2008, respectively. 

'https : //safe.nrao . edu/wiki/bin/view/CICADA/NGNPP 
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Figure 2 Optical spectrum for J 1023, taken on Dec. 23 and 25, 2008. Flux was normalized by 
dividing by the best-fit continuum spectrum. Note the absence of emission lines; this spectrum 
is consistent with that observed in quiescence and not consistent with that observed during 
J1023's active phase. 
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The pulse profile next to (C) is based on a 1400-MHz band observation taken at the Arecibo 
Observatory on Nov. 18, 2008 and coherently dedispersed with the ASP backend (33), while 
the others are based on the same observations as the grayscale plots and generally suffer from 
poorer time resolution. Smearing in the profiles obtained by coherent dedispersion is due to 
pulse drift and inaccurate DM values; the time and frequency subintegrations have been aligned 
with an improved ephemeris but smearing within each subintegration remains. 



Obtaining a timing solution: From each observation of PSR J 1023+0038, we produced 
folded profiles for time and frequency subintegrations. We then established a reference profile 
for each combination of band and backend, usually from the best observation available with that 
configuration, though we re-used the standard profile from the 1400 MHz ASP observations for 
the 1388 MHz and 1788 MHz GASP observations. We produced a pulse time-of-arrival (TOA) 
for each subintegration (generally 30 s or 60 s) by cross-correlating the reference profile with 
the folded profile for that subintegration §34\ . These arrival times are based on observatory time 
standards which were corrected to UTC using the GPS satellites. We fed these TOAs into the 
pulsar timing program TEMPC^to obtain a timing solution, allowing for an arbitrary time offset 
of each band/backend combination. The solution we obtain is referenced to the DE405 solar 
system ephemeris. 

The initial timing solution we obtained showed the timing effects of the DM variations we 
observed, as well as possibly some other systematics. We excluded any data that showed signs 
of short-term DM variations, including all data between orbital phases 0.15 and 0.55, as well 
as most low-frequency data. Specifically, the only data below 1 GHz we used was one pair of 
800 MHz GBT observations separated by a day and one 350 MHz GBT observation. These 
data sets did not show signs of significant short-term DM variations, and the inter-instrument 
jump we assumed absorbed any effects of long-term DM variations. To mitigate any remaining 
systematics, we artificially increased the uncertainties on each individual data set until it was fit 
with a reduced x 2 of approximately one. 

We have attached the resulting lists of pulse arrival times in TEMPO format to this support- 
ing online material. 



Estimating DM variations: The global DM value we report is obtained from the Oct. 18, 
2008 GUPPI observation. We produced TOAs for four frequency subbands for each time subin- 
tegration, aligning the folded data with a common template, and used the pulsar timing package 
tempo to fit for the DM that minimized the weighted residuals. This process yields a very 
precise DM, but when we applied an analogous procedure to several other data sets, the results 
varied from 14.32 to 14.33 pc cm~ 3 . We applied a similar process to the data taken at Arecibo 
with the WAPPs, obtaining a DM value for each time subintegration. These values showed 
the large DM excesses around orbital phase 0.3, but they also confirmed the variability of the 

"http : / /www .atnf.csiro.au/ research/pulsar/ tempo/ 
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average DM, giving values stable within an orbit but varying from orbit to orbit by as much as 
0.01 pc cm -3 . 
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Table 2. 



Published observations of J 1023. 



Date Type Reference Comment 



1 QOT Mra; 7 

lyyZ. INOV / 


436 MHz pulsar survey 3 I 


35 


i No detection (< 3 mjy) 


1770 ^VUg J 


1 .4 GHz continuum b 




< 1.8 mJy 


1 QQS Ann fi 

iyyo /vug o 


1 .4 GHz continuum b 


8 


< 3.4 mJy 


1 QQS A hit in 


1 .4 GHz continuum 13 


8 


6.56 mJy 


1 QQQ Mnr 01 
lyyy IVldJ 


u', g', r', i', z' photometry 1 


5o| G star color 


ZUUU May o 


Spectroscopy' 1 


S 


Blue, emission lines 


zuuu iviay y 


Spectroscopy ESO archive Blue, emission lines 


ZUUU iNOV Zl 


Fast photometry* 


8 


Flickering by ~ 1 mag 


onnn xt™/ to 

ZUUU iNOV ZZ 


Fast photometry* 


8 


Flickering by ~ 1 mag 


zuuu uec zz 


Fast photometry* 


8 


Flickering by ~ 1 mag 


zuuu Dec Zj 


Fast photometry* 


8 


Flickering by ~ 1 mag 


2000 Dec 24 


Fast photometry* 


8 


Flickering by ~ 1 mag 


2001 Feb 1 


Spectroscopy I 


36 


i Blue, emission lines 


2001 Uec 1U 


Spectroscopy 5 


9 


Blue, double-peaked emission lines 


2002 May 1 1 


Spectropolarimetry* 1 


21 


G star spectrum, no significant polarization 


2002 Dec 30 


Fast photometry 1 


27 


Smooth, low-amplitude flickering 


2003 Oct 16 


8.4 GHz continuum 13 


37 


No detection 


2003 Jan 28 


Fast photometry* 


27 


Smooth low-amplitude flickering 


2003 Jan 29 


Prist nhntnmptreJ 


27 


Smfiiiln Inw-nmnlitnHp. flickprmo 

LJlllvJUlll, 1UW LllllUllLLlLlC lIlLlvCllllw 


2003 Jan 30 


Pnst ntintnmptrvJ 


27 


Smfiiiln Inw-nmnlitiiHf flirlfpnno 

LJlllvJUlll, 1UW LllllUllLLlLlC lIlLlvCllllw 


2003 Jan 31 


Fast photometry' 


27 


Smooth, low-amplitude flickering 


2003 Feb 03 


Fast photometry' 


27 


Smooth, low-amplitude flickering 


2003 Feb 05 


Fast photometry! 


27 


Smooth, low-amplitude flickering 


2003 Feb 21 


Fast photometry' 


27 


Smooth, low-amplitude flickering 


2003 Feb 24 


Fast photometry' 


27 


Smooth, low-amplitude flickering 


2004 Feb 16 


Spectropolarimetry 15 


21 


G star spectrum, no significant polarization 


2004 May 12 


X-ray 1 


21 


Hard power-law spectrum 


2004 May 12 


X-ray m 


21 


Irregular, 60% modulation 


2004 May 12 


B photometry 11 


21 


Smooth periodic 40% modulation 


2004 May 13 


Photometry 


21 


Smooth periodic 30% modulation 


2004 May 23 


Spectroscopy 13 


21 




2003 Jan 31 


Spectroscopy* 3 


7 




2004 Jan 18 


Spectroscopy* 3 


7 




2004 Jan 19 


Spectroscopy* 3 


7 




2004 Jan 20 


Spectroscopy* 3 


/ 




2004 Mar 8 


Spectroscopy* 3 


7 




2004 Mar 9 


Spectroscopy* 3 


7 




2004 Nov 18 


Spectroscopy* 3 


7 




2004 Feb 29 


Filter photometry* 3 


7 




2004 Feb 29 


Fast / photometry* 3 


7 




2004 Mar 1 


Fast / photometry* 3 


7 




2004 May 12 


Fast B, V, I photometry q 


/ 




2004 May 16 


Fast B, V, I photometry 01 


7 




2004 May 20 


Fast B, V, I photometry ! 


7 






H a , I' , r' imaging 1 I 


18 


i No H a excess 



Note. — Horizontal lines indicate the beginning and end of the active phase. 
a The Australia Telescope National Facility Parkes radio telescope 
b The Very Large Array, operated by the National Radio Astronomy Observatory 
°The SDSS 2.5 m telescope at the Apache Point Observatory 
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d Lick Observatory 3 m telescope 

°ESO Faint Object Spectrograph and Camera on the 2.2 m telescope at La Silla 
f Kitt Peak National Observatory 2.1 m 

g Double-imaging spectrograph on the Apache Point Observatory 3.5 m telescope 
h Kitt Peak Bok Telescope, the spectropolarimeter SPOL 

'The 1.9 m telescope at the Sutherland site of the South African Astronomical Observatory 
j The 1.0 m telescope at the Sutherland site of the South African Astronomical Observatory 
k The 6.5 m Multi-Mirror Telescope 
'XMM-Newton EPIC-pn 
m XMM-Newton EPIC-MOS1/2 
n XMM-Newton OM 

°US Naval Observatory Flagstaff Station 1 m telescope 
p The Hiltner 2.4 m telescope 
q The 1.3 m McGraw-Hill telescope 
r The 2.5 m Isaac Newton Telescope 
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Table 3. Our Observations of J1023. 



Date 


MJD range 


Orbital phase 


Instrument 


Frequency (MHz) 


Comment 


2007 Jun 25 


54276.00 


0.85 


GBT Spigot a 


350 


Nearby track 


2007 Jun 28 


54279.98 


0.96 


GBT Spigot 3 


350 


Discovery 


2008 Oct 18 


54757.54 


0.71 


GBT Spigot*-" 


350 




2008 Oct 18 


54757.58-54757.66 


0.91-1.32 


GBT GUPPI* /Spigot b 


2000 




2008 Oct 18 


54757.68 


0.40 


GBT GUPPI/Spigot b 


2000 




2008 Oct 18 


54757.74 


0.69 


GBT GUPPI/Spigot b 


2000 




2008 Oct 21 


54760.76-54760.79 


0.95-1.09 


GBT Spigot*-" 


800 




inns rw io 
ZUUo wet Zz 


J4 / Ol .4 /— J4/01 . J 1 


n ^7 n 77 

U.J /— u. / j 


uis i spigot ' 


snn 
5UU 




2008 Oct 22 


54761.53-54761.55 


0.83-0.96 


AO ASP/WAPPs c 


327 




2008 Oct 24 


54763.82-54764.06 


0.39-1.59 


Parkes 10 cm* /50 cm d 


3000/700 




2008 Oct 27 


54766.53-54766.54 


0.09-0.14 


AO ASP* AVAPPs c 


1400 




2008 Nov 1 


54771.50-54771.68 


0.20-1.10 


GBT Spigot/GUPPI/GASP* 


<f 1400 




2008 Nov 6 


54776.49-54776.50 


0.37-0.43 


AO ASP/WAPPs c 


1400 




2008 Nov 1 1 


54781.03-54781.25 


0.30-1.39 


WSRT PuMa-HS 


150 




2008 Nov 12 


54782.47-54782.50 


0.55-0.70 


AO ASP* AVAPPs e 


1400 




2008 Nov 18 


54788.45-54788.48 


0.77-0.91 


AO ASP* AVAPPs e 


1400 




2008 Nov 24 


54794.43-54794.46 


0.92-1.08 


AO ASP* AVAPPs c 


1400 




2008 Dec 2 


54802.42-54802.44 


0.28-0.40 


AO ASP/WAPPs 6 


1400 




2008 Dec 12 


54812.39-54812.42 


0.62-0.74 


AO ASP* AVAPPs e 


1400 




2008 Dec 22 


54822.36-54822.39 


0.95-1.09 


AO WAPPs h 


1400 




2008 Dec 31 


54831.09-54831.33 


0.03-1.22 


WSRT PuMa-II 1 


350 




2009 Jan 2 


54833.33-54833.36 


0.31-0.48 


AO ASP/WAPPs 


1400 




2009 Jan 5 


54836.16-54836.33 


0.60-1.48 


GBT GASP* /GUPPP 


1700 




2009 Jan 16 


54847.28-54847.30 


0.76-0.84 


AO ASP* AVAPPs c 


1400 





*Used in the timing solution. 

a We used the Spigot with a 50 MHz band with 2048 frequency channels and a time resolution of 81.92 /is. 

b We used the Spigot with an 800 MHz band with 2048 frequency channels and a time resolution of 81.92 /is. GUPPI recorded full Stokes 
parameters for an 800 MHz band with 2048 frequency channels and a time resolution of 40.96 fj.s. In both cases only about 600 MHz of the 
bandwidth is usable. 

c We used two WAPPs, each with a 12.5 MHz band with 512 frequency channels and a time resolution of 64 /is. ASP recorded a bandwidth 
of 16 MHz centered at 327 MHz. 

d The dual-band receiver synthesized two filterbanks, one 64 MHz wide with 256 channels and centered at 590 MHz, and the other 768 MHz 
wide with 256 channels centered at 3032 MHz. Both had a time resolution of 80 fj.s. 

e We used three WAPPs, with center frequencies 1170, 1370, and 1570 MHz, using bandwidths of 100 MHz with 256 channels and a time 
resolution of 64 fxs. ASP recorded a 52 or 56 MHz bandpass centered on 1412 MHz. 

'We used the Spigot with an 800 MHz band with 2048 frequency channels and a time resolution of 81.92 /is. GUPPI recorded full Stokes 
parameters for an 800 MHz band with 2048 frequency channels and a time resolution of 40.96 /zs. In both cases only about 600 MHz of the 
bandwidth is usable. GASP recorded an 84 MHz bandpass centered at 1388 MHz. 

g We used the PuMa-II backend to coherently dedisperse 8 2.5 MHz channels. 

h We used the same three-WAPP configuration as in f above, but the ASP data was accidentally corrupted. 
'We used the PuMa-II backend to coherently dedisperse 6 10 MHz channels. 

j We used GUPPI in an online folding mode, recording full Stokes parameters for an 800 MHz band with 2048 frequency channels and a 
time resolution of 40.96 /is; only 600 MHz of the bandwidth was usable. GASP recorded a 76 MHz bandpass centered at 1788 MHz. 
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